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We point out that the extra neutralinos and charginos generically appearing in a large class of 
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from the Minimal Supersymmetric Standard Model (MSSM). These signatures involve longer decay 
chains, more visible transverse energy, higher multiplicities of jets and leptons, and less missing 
transverse energy than in the MSSM. These features make the gluino harder to discover for certain 
types of conventional analyses, for example in the 6 jet channel we show that a 1 TeV MSSM 
gluino may resemble an 800 GeV Ee gluino. However, the enriched lepton multiplicity provides 
enhanced 3- and 4-lepton signatures, which are much more suppressed in the case of the MSSM. As 
an example we analyse various scenarios in the EgSSM, demonstrating the utility of a CalcHEP model 
that we have now made publicly available on HEPMDB (High Energy Physics Models DataBase). 
After extensive scans over the parameter space, we focus on representative benchmark points. We 
perform detailed Monte Carlo analyses, concentrating on the 3-lepton channel at the 7, 8, and 14 TeV 
Large Hadron Collider (LHC), and demonstrate that inspired models are clearly distinguishable 
from the MSSM in gluino cascade decays. We emphasise to the LHC experimental groups that 
the distinctive features present in gluino cascade decays, such as those discussed here, not only 
represents a unique footprint of a particular model but also may provide the key to an earlier 
discovery of supersymmetry. 
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I. INTRODUCTION 

The general idea of softly broken supersymmetry (SUSY) provides a very attractive framework for physics models 
beyond the standard model (BSM), in which the hierarchy problem is solved and the unification of gauge couplings 
can be realised [JJ. The simplest such extension, the Minimal Supersymmetric Standard Model (MSSM), has been 
subjected to particularly close scrutiny at the CERN Large Hadron Collider (LHC). For example, searches at ATLAS 
Ej and CMS OH], under certain assumptions, constrain the MSSM gluino mass to be greater than about 500-900 
GeV. This mass limit range reflects the fact that the precise limit depends crucially on the assumptions one makes 
about the rest of the SUSY spectrum, in particular the mass spectrum of the squarks, charginos, and neutralinos. As 
we emphasise in this paper, mass limits and the LHC discovery potential also depend on the particular SUSY model 
under consideration and the quoted mass limits strictly only apply to the constrained MSSM or simplified models of 
SUSY[5] with very few particles. In other, well motivated, non-minimal SUSY extensions of the SM, such as those 
discussed below, the enriched particle content alters the SUSY discovery potential at the LHC and this should be 
taken into account in the phenomenological and experimental studies. 

Despite its many attractive features, the MSSM suffers from the \i problem. The superpotential of the MSSM 
contains the bilinear term /iHdH u , where Hd and H u are the Higgs doublet superfields. In order to get the correct 
pattern of electroweak (EW) symmetry breaking the parameter fi is required to be in the TeV region. At the same 
time, the incorporation of the MSSM into grand unified theories (GUTs) or string theory implies that /i could be of the 
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order of the GUT or Planck scales, or possibly zero. None of these possibilities is phenomenologically and theoretically 
acceptable. The Next-to-Minimal Supersymmetric Standard Model (NMSSM) [5] attempts to address this problem by 
postulating an additional gauge singlet superfield 5 with the interaction XSHdH u in the superpotential together with 
an 5 3 interaction in order to break an accidental global U(l) symmetry to a discrete Z 3 symmetry. At low energies 
(~ TeV) the scalar component of 5 acquires a non-zero vacuum expectation value (VEV) s, giving an effective fj, term. 
However the resulting Zj, discrete symmetry is broken at the same time, leading to potentially dangerous cosmological 
domain walls [6] . 

A solution to the /z problem which does not suffer from the problems of the NMSSM arises within Eq inspired 
SUSY models. At high energies the Eq gauge symmetry arising from GUTs or string theory in extra dimensions can 
be broken to the rank-5 subgroup SU(3) C x SU(2)l x U(1)y X 1/(1)', where in general 

U(l)' = U(l) x cos 6 + U(l)j, sin0 (1.1) 

and the two anomaly-free 17(1)^, and U(l) x symmetries originate from the breakings Eg —> 50(10) x U(l)^, then 
50(10) — > SU(5)x U(l) x . The extra U(l)' gauge symmetry forbids 5 3 interactions since the SM-singlet superfield 
5 is now charged under U(l)'. In addition the bilinear /i term is also forbidden if 6 ^ or ir, while the interaction 
XSHdH u is allowed in the superpotential. 

At low energies (~ TeV) the scalar component of 5 acquires a non-zero VEV s breaking U(l)' and giving rise to a 
massive Z' gauge boson together with an effective fj, term. Within the class of rank-5, Eq inspired SUSY models with 
an extra 17(1)' gauge symmetry, there is a unique choice of the extra Abelian gauge group that allows right-handed 
neutrinos to be uncharged. This is the U(1)n gauge symmetry given by 9 = arctan vl5 [71 [5]. In this particular case, 
called the Exceptional Supersymmetric Standard Model (E 6 SSM), based on the SU(3) C x SU(2) L x U(1)y x U(1)n 
gauge group, the right-handed neutrinos may be superheavy, allowing a high scale see-saw mechanism in the lepton 
sector [7J|H]. 

An elegant feature of the Eq inspired models with an extra Z' gauge boson at the TeV scale is that the conditions 
of anomaly cancellation may be satisfied by matter content filling out complete 27 representations of Eq surviving to 
the TeV scale. One can therefore have three 27 representations, with each containing a generation of SM matter and 
fields with the quantum numbers of Higgs doublets and SM-singlets. Thus such models predict, in addition to three 
SUSY families of quarks and leptons, also three SUSY families of Higgs doublets of type H u , three SUSY families of 
Higgs doublets of type Hd, three SUSY SM singlets Si, and three SUSY families of charge ^fl/3 colour triplet and 
anti-triplet states D and D, which can get large mass terms due to the (third) singlet VEV (53) = s/y/2 which is also 
responsible for the \x term. 

In this paper we study gluino production and decay within the E 6 SSM, expected to be the main discovery channel 
for this model. In fact the analysis also applies to a larger class of Eq models in which the matter content of three 27 
representations of Eq survives to the TeV scale. We study the EgSSM as a concrete example and at the same time 
demonstrate the use of the E6SSM-12.02 CalcHEP model that we have now made publicly available on HEPMDB. 
After extensive scans over the parameter space we create a number of representative benchmarks for different E 6 SSM 
scenarios, considering Z' and Higgs boson physics, the LSP dark matter relic density and direct detection cross-section, 
and the perturbativity of dimensionless couplings. We then analyse representative MSSM and EgSSM benchmarks 
consistent with the particle recently discovered at the LHC [9l [10] being the SM-like lightest Higgs boson. 

The analysis is motivated by the fact that the gluinos are the expected to be the lightest strongly interacting particles 
in Eq models |llj . so one should expect them to have the largest production rate. In particular, we are interested 
in gluino cascade decays which not only provide observable signatures, but are also important for distinguishing the 
Eq inspired models from the MSSM. As we shall see, there are important differences between the two cases, which 
can affect the respective search strategies. The differences arise because of the extra Higgsino and singlino states 
predicted to be part of the Eq matter content, as above. In particular, in the Eq inspired models, relative to the 
MSSM, there are two extra families of Higgsinos, H" and H^, together with two extra singlinos, S a , where a = 1,2^] 
There is also a third singlino, 5 3 , similar to the NMSSM singlino, which mixes with the bino', both states having a 
large mass of order the effective fi parameter. The remaining extra Higgsinos and singlinos may be lighter than the 
gluino. Indeed, it is possible to show that at least two linear combinations of the states H", H^, S a must be lighter 
than or of order Mz/2. If these states mix with the usual neutralinos of the MSSM or NMSSM then the lightest 
supersymmetric particle (LSP) will inevitably be one of these states, leading to longer decay chains. For example, in 
regions of MSSM parameter space where the bino is the LSP the gluino typically undergoes a cascade decay to the 
bino. In the Eg inspired models the bino will mix with the extra Higgsinos and singlinos and the predominantly bino 



1 Note that the first and second family of Higgs doublet and singlet fields H", H^, and S a predicted by the EgSSM do not develop VEVs 
and are called "inert". 
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state will subsequently decay into some lighter state having a mass of order Mz/2, thereby typically giving a longer 
gluino cascade decay chain and producing less missing energy due to the lighter mass of the LSP0 For simplicity we 
shall assume that the D and D states, as well as the NMSSM type singlino S 3 , are all heavier than the gluino and 
so are irrelevant for gluino cascade decays. Similarly we shall also assume all squarks and sleptons and Higgs scalars 
and pseudoscalars (with the exception of the SM-like Higgs boson) to be heavier than the gluino. These assumptions 
are motivated by the parameter space of the constrained E§ inspired models |llH13j . 

The main result of our analysis is that in Eq inspired models the gluino decays into the LSP with longer chains 
which involve more jets and leptons and less missing energy than in the MSSM. This happens because the would-be 
LSP (e.g. an MSSM-like bino dominated neutralino) undergoes further decays to the extra light neutralinos and 
charginos predicted by the E§ inspired models. As a result, the characteristics of the signal, such as lepton and jet 
multiplicity, missing transverse momentum, effective mass, etc., are altered in the EgSSM as compared to the MSSM 
case even after the matching of the gaugino masses in both models. Therefore, the search strategies designed for the 
MSSM need to be modified for the EgSSM case, while one should stress that the gluino mass limits for the MSSM 
are not applicable to the E 6 SSM gluinos. 

The layout of the remainder of this paper is as follows: In section II we summarise the relevant features of the 
EgSSM, its theoretical constraints and its experimental constraints not relating to gluino detection. We also present 
the results of parameter scans and discuss the viable parameter space, before going on to introduce a set of benchmark 
points, including those which form the basis of our analysis. The model implementation is discussed in section III 
(with a more in-depth description of the publicly available code in Appendix B). We then present our analysis of the 
different signatures and search prospects for different strategies for each model in section IV. The conclusions are in 
section V. 



II. MODEL SETUP AND PARAMETER SPACE 



Before going on to consider the prospects for the production and detection of gluinos, which are the main focus of 
this paper, discussed in section IV, we must first determine the limits on the EgSSM from other experimental, and 
cosmological, considerations. We begin by giving an introduction to the EgSSM (subsections II A| and II B I, explaining 
the important features. We then discuss experimental constraints relating to Z' and Higgs boson physics and to exotic 
coloured particles (subsection II C ) , before going on to discuss dark matter considerations (subsection II D[ ) . In light 
of these discussions we then show the results of some parameter space scans (subsection HE) and produce a set 
of benchmark points (subsection II F). These benchmarks example various viable scenarios for the EgSSM and we 
explain their features and issues. Although the benchmarks presented look very different from each other from many 
points of view, it turns out that they look very similar in terms of their gluino decay signatures. In section IV we 
therefore mostly show results for two particular benchmarks that demonstrate the qualitative differences between the 
MSSM and E G models for gluino searches. 



A. The E 6 SSM 



As discussed in the Introduction, the EgSSM involves a unique choice for the extra Abelian gauge group that allows 
zero charges for right-handed neutrinos, namely U(l)jf. To ensure anomaly cancellation the particle content of the 
EgSSM at the TeV scale is extended to include three complete fundamental 27 representations of Eq, apart from the 
three right-handed neutrinos which are singlets and do not contribute to anomalies and so may be much heavier. The 
27 representations of E§ decompose under the SU(5) x U(1)n subgroup of E6 as follows: 

27, -> (10, 1), + (5* , 2), + (5*, -3)< + (5, -2), + (1, 5), + (1, 0), . (ILL) 

The first and second quantities in brackets are the SU(5) representation and extra U(1)n charge (where a GUT 



normalisation factor of l/v40 is omitted OH]) respectively, while i is a family index that runs from 1 to 3. An 
ordinary SM family, which contains the doublets of left-handed quarks Qi and leptons Lj, right-handed up- and 
down-quarks («? and d|), and right-handed charged leptons is assigned to (10, l)j + (5*, 2)j. Right-handed neutrinos 



Nf should be associated with the last term in Eq. (II. 1), (1,0)^ The next-to-last term, (1,5)*, represents SM-singlct 
fields Si, which carry non-zero U(1)n charges and therefore survive down to the EW scale. The pair of SU(2)l- 
doublets {Hf and Hf) that are contained in (5*, —3), and (5,— 2)j have the quantum numbers of Higgs doublets. 



2 The decay to bino is expected to happen before the subsequent decay into a lighter state since these lighter states are expected to have 
small mixing to the MSSM-like sector, for reasons explained in section II. 
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They form either Higgs or inert Higgs SU(2)l multiplets. Other components of these fundamental SU(5) multiplets 
form colour triplet and anti-triplet of exotic quarks Di and Di with electric charges —1/3 and +1/3 respectively. 
These exotic quark states carry a B — L charge +2/3, twice that of ordinary quarks. In addition to the complete 
27^ multiplets the low energy matter content of the E 6 SSM can be supplemented by an SU{2)i 1 doublet L4 and 
anti-doublet L4 from extra 27' and 27' representations to preserve gauge coupling unification [14 . These states will 
typically be much heavier than the gluino and so will play no role in the present analysis, although they may play a 
role in leptogenesis [15] . 

As in the MSSM the gauge symmetry in the EgSSM does not forbid lepton and baryon number violating operators 
that result in rapid proton decay, although the situation is somewhat different. The renormalisable superpotential of 
the EgSSM automatically preserves i?-parity because of the enlarged gauge symmetry, however, although the B — L 
violating operators of the i?-parity violating MSSM are forbidden by the gauge symmetry, there are new B and L 
violating interactions involving the new exotic particles. In the EgSSM these terms are removed by imposing an 
exact Z2 symmetry on the superpotential. There are two options for this symmetry under which the exotic quarks 
are either interpreted as diquarks or leptoquarks, leading to B and L conservation. Furthermore, the extra particles 
present in Eg inspired SUSY models give rise to new Yukawa interactions that in general induce unacceptably large 
non-diagonal flavour transitions. To suppress these effects in the EgSSM an additional approximate Z^ symmetry is 
imposed. Under this symmetry all superfields except one pair of H d and B % u (say Hd = i?5 and H u = H 3 ) and one 
SM-type singlet field (S = 5* 3 ) are odd. Ignoring L4 and L4, the Z^ symmetry reduces the structure of the Yukawa 
interactions to 

We 6 ssm — ^S(H u Hd) + \ a pS(H2Hu) + fuapS a {H d H u ) + fdafjS a {HdH^) 
+ hYj(H u Qi)Uj + h^HdQi)^ + h%(HdU)e] + h%(H u Li)N? 

+ ^MijftfNj + KiSibibi) , (II.2) 

where hats denote superfields, a, f3 — 1,2 and i, j = 1, 2, 3. The inert Higgs doublets and SM-singlets have suppressed 
couplings to matter, whereas the third generation SU(2)l doublets H u and Hd and SM-type singlet field S, that are 
even under the Z% symmetry, play the role of Higgs fields, radiatovely acquiring VEVs. At the physical vacuum their 
scalar components develop VEVs 

W vi \ ,„ s 1 ( A lrn s 



W-^SJ. W-^UJ' <S> = 7!' (IL3) 

generating the masses of the quarks and leptons. Instead of v\ and V2 it is more convenient to use tan/3 = i^/vi and 
v — \Jv\ + v\ — 246 GeV. The VEV of the SM-singlet field, s, breaks the extra U(1)n symmetry , generating exotic 
fermion masses and also inducing a mass for the Z' boson. Therefore the singlet field S must acquire a large VEV in 
order to avoid conflict with direct particle searches at present and past accelerators. 

We shall assume that the Kj are sufficiently large that the exotic Di,Di states are much heavier than the gluino 
and so will play no role in gluino decays. Similarly the right-handed neutrinos Ni will be neglected since they are 
assumed to be very heavy. 



In this paper we shall be mainly concerned with the first line of Eq. (II.2), namely the couplings of the form 
S % HiHl^ > where k — 1,2,3 label the three families of Higgs doublet and singlet superfields predicted in the 
EgSSM. In particular we shall be concerned with the resulting chargino and neutralino mass terms coming from 
such couplings involving one third-family scalar component and two fermion components, i.e. SH^H^, S l HdH^, and 
S l H° d H u . In the EgSSM the presence of the extra Higgsinos and singlinos H", H2, and S a means that the chargino 
and neutralino mass matrices are extended, as discussed in the following subsection. 



B. Neutralino and chargino mass matrices 

In the MSSM T] there are four neutralino interaction states, the neutral wino, the bino, and the two neutral 
Hig gsinos. In the USSM |16j . a model similar to the NMSSM, but where the t^(l) introduced by adding the singlet 
field S is gauged instead of reduced to a Z3, two extra states are added, the singlino and the bino'. In the conventional 
USSM interaction basis we write the six neutralinos as the column vector 

XL = (B W 3 H° U \S B' f. (II.4) 
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Neglecting bino-bino' mixing (assumed to be small as justified in Ref. 
basis becomes 



the USSM neutralino mass matrix in this 



MjSSM 



Mi 





—mzswcp 


mzSy/Sp 










-mzswcp 


M 2 
m z c w Cf3 


mzcwcp 



-m z c w sp 
-V 








m z swSj3 


-mzcwsp 


-V 















9'iVC(sQ% 






g'isQ? 


g'isQ 1 : 

M[ 



(11.5) 



where M±, M2, and M{ are the soft gaugino masses and fx s — \v/\2. The variables s^ny* and Cp(y/) stand for sin 
and cos of /3(the Weinberg angle) and Q^ d u ™ are the (H^, H u , S) U(1)n charges (—3, —2, 5)/v40. 

In the EgSSM the neutralino sector is extended to include the additional six neutral components of H™, H2, and 
S a , where a = 1,2. We then take the full list of twelve neutralino interaction states into account in column vector 



Xint 



B W 3 H% H° I S B' I H% 2 H° 2 S 2 I H^i H°i &i ) T 



(II.6) 



The first four states are the MSSM interaction states, the S and B' are the extra states added in the USSM, and the 
final six states are the extra inert doublet Higgsinos and Higgs singlinos that come with the full E 6 SSM. Under the 
assumption that only the third generation Higgs doublets and singlet acquire VEVs the full Majorana mass matrix is 
then [T7] 



M n 

JW E 6 SSM 




(II.7) 



where the sub-matrices involving only the inert interaction states are given by 



A, 



a/3 



A Q(3 s fu/3 a v sin /3 

Xp a S fdfSaVCOS/3 

f ua fSV Sill (3 f dal 3V COS /3 



(II. 



and the Z 2 breaking sub-matrices are given by 



B n 



1 

71 



/ 



Xd a s z a vsinj3 

v cos /3 

x ua v sin j3 XdaV cos /3 

V 



(II.9) 



and involve the small Z 2 violating Yukawa couplings that were neglected in Eq. ( II. 2 ) , x ua , Xda, and z a . Since these 
coupling are small, 

the inert neutralino sector is only weakly coupled to the USSM sector, and may be considered separately from it 
to good approximation. However we emphasise that the Z 2 violating couplings are essential in order for the lightest 
neutralino from the USSM sector to be able to decay into inert neutralinos and that these couplings are not expected 
to be zero. Exact Z 2 would also render exotic D and D states stable. 

Similarly, we take our basis of chargino interaction states to be 



Xint 
Xint 



where 



Xint 



and 



H 



,12 



(11.10) 
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The corresponding mass matrix is then 



A I, 



E 6 SSM 



c 



where 



C = 



( Mi V2m w sin /? 
\2mw cos j3 fi 











1 Xd2S \x d \8 



V 



X U 2S 
X U \S 



J 2 



\22S 



71 A2lS 



(11.11) 



It is clear that a generic feature of the EgSSM is that the LSP is usually (naturally) composed mainly of inert 
singlino and ends up being typically very light. One can see this by inspecting the new sector blocks of the extended 
neutralino mass matrix in Eq. (II.7I, such as An, and assuming a hierarchy of the form \ a ps 3> fuapv, fdupv- This 



is a natural assumption since we already require that s 3> v in order to satisfy the current experimental limit on the 
Z' mass of around 2 TeV, as discussed below and for example in Ref. [T5] . 

We emphasise again that for both the neutralinos and the charginos we see that if the Z^ breaking couplings are 
exactly zero then the new part of the EgSSM neutralino mass matrix becomes decoupled from the USSM mass matrix. 
However, although approximate decoupling is expected, exact decoupling is not, and will therefore not be considered. 



C. Experimental constraints 

The most recent limit on the U(1)n Z' mass, set by the CMS [T^j, searching for dilepton resonances, is mz' > 
1800 GeV at a confidence level of 95%. Although the limit on the mass of the Z' boson associated with the extra 
U(1)n of the EgSSM can be inferred from this analysis, this analysis neglects any other matter beyond that of the 
SM. When decays of the Z' boson into inert neutralinos (inert Higgsino and singlino dominated mass eigenstates) 
are considered the Z' width tends to increase by a factor of about 2 (see for example Ref. [13] , although we confirm 
the result in our analysis). This then means that the branching ratio into leptons is decreased by a factor of about 
2. Estimating the effect of halving this expected branching ratio on the analysis in Ref. [H] one can read off a 95% 
confidence level lower bound of around 1600 GeV. This implies s > 4400 GeV. 

When the singlet VEV is this large, as required, the Higgs boson spectrum typically becomes rather hierarchical 
with a lightest mass eigenstate that participates in interactions as a SM-like Higgs boson and much heavier Higgs 
boson states that are approximately decoupled. Indeed the Higgs candidate recently discovered at the LHC can easily 
be accommodated in the EgSSM [20] ■ Although in the MSSM large loop corrections are required in order for this 
limit to be satisfied, in the EgSSM this is easier to achieve since there are extra contributions to the SM-like Higgs 
mass at tree-level due to U(1)n D-terms. 

If the exotic diquarks (or leptoquarks) are light enough 

they would produce spectacular signatures at the LHC and these exotic states already have strong limits on their 
masses. The E e diquarks are excluded for masses below 3.5 TeV [21] • Since these particles' masses must be so large 
they do not play a role in gluino cascade decays and are excluded from our analysis. 



D. Dark matter considerations 



Stringent constraints on the EgSSM inert parameter space come from considerations relating to dark matter. In 
the EgSSM as described thus far the LSP is typically one of the two necessarily light states from the inert neutralino 
sector. As such this inert neutralino LSP becomes a dark matter candidate. The EgSSM has been previously studied 
as a model attempting to explain the observed amount of thermal relic cold dark matter [T7[ [22] . Unfortunately this 
dark matter scenario is now severely challenged by the most recent XENON100 dark matter direct detection limits 
[23j . The reason is essentially as outlined in the following paragraph and a more detailed analysis can be found in 
Ref. [22]. 

In the EgSSM the LSP is generically singlino dominated, a situation which arises from the extended neutralino 
mass matrix in Eq. (II. 7 1 under the condition s ^> v. One can show that [T7[ [22] if there is no hierarchy in the 
Yukawa couplings then the LSP would not annihilate very efficiently at the time of thermal freeze-out and would 
therefore lead to an unacceptable overabundance of dark matter in the Universe. On the other hand, by allowing the 
largest fda/3 and f U a/3 couplings to be significantly larger than the largest X a p couplings and tan fi to be less than 
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about 2 the observed amount of dark matter can be predicted. In this case the LSP is heavier and, although still 
inert singlino dominated, has substantial inert Higgsino admixtures and can annihilate efficiently enough in the early 
universe via an s-channel Z boson. The largest \ a p coupling cannot be too small, otherwise the inert charginos would 
be too light to have so far escaped detection. At the same time the largest fdap and f U a/3 couplings cannot be too 
large if it is required that perturbation theory remains valid up to the GUT scale. This being the case, the LSP and 
the NLSP cannot be made much larger than about 60 GeV [55]. In this dark matter scenario there should be some 
suppression of the coupling of the LSP to the Z boson, by partial cancellation between the up-type and down-type 
inert Higgsino components, in order to be consistent with the precision measurement of the invisible Z boson width 
from LEP, if the LSP mass is below half of the Z mass. For given Yukawa couplings increasing tan/3 has the effect 
of both suppressing the LSP mass and increasing the coupling to the Z boson, lessening this cancellation. Although 
this partial cancellation can occur in the coupling of the LSP to the Z boson, the coupling of the LSP to the SM-like 
Higgs boson is necessarily large if the LSP is to produce the observed amount of dark matter (due to its inert Higgsino 
admixtures) . This in turn means a large spin-independent direct detection cross-section, larger than is now consistent 
with experiment. However, if the relic abundance is less than the observed value, then the direct detection constraint 
can be avoided, and that is the strategy that we follow for benchmark points of this kind, as we now discuss in more 
detail. 

It is important to note that for both the MSSM and the standard EgSSM we analyse points where less than 
the observed amount of dark matter is predicted and assume that the majority of dark matter is not made up of 
MSSM/EgSSM neutralinos. This choice of the parameter space is actually dictated by limits from direct detection 
experiments: if less than the total amount of the dark matter in the universe is made up of LSPs, then the expected 
number of events for a direct detection experiment for a given LSP direct detection cross-section will be correspondingly 
smaller. In the EgSSM such points require the previously discussed hierarchy of Yukawa couplings appearing in the 
inert block of the neutralino mass matrix and in this case if the LSP mass is below half of the Z mass then tan/3 
still cannot be too large in order for the LSP not to contribute too much to the invisible width of the Z boson, as 
measured at LEP. Large values of tan j3 mean that the up-type inert Higgsino admixture in the LSP greatly outweighs 
the down-type inert Higgsino admixture, necessarily leading to a too large coupling to the Z boson. 

Alternatively, in a variation of the model known as the EZSSM [21] the situation is quite different and the model 
may be responsible for all dark matter and consistent with all experiment. Here an extra discrete symmetry is 
imposed that forbids the terms in the superpotential involving the inert singlet superfields, i.e. the f ua p and fd a p and 
Z2 violating x ua and Xda Yukawa couplings are forced to be zero. This then means that the inert singlinos are exactly 
massless and decoupled from the rest of the neutralinos. In this variation of the model the lightest non-inert-singlino 
LSP is absolutely stable. 

If this stable particle is the bino and there are a pair of inert Higgsinos close by in mass then the bino can be 
responsible for all of dark matter [24] . The massless inert singlinos themselves slightly increase the expansion rate of 
the universe prior to nucleosysnthesis in agreement with observation of the 4 He relic abundance and in this scenario 
there would currently be a cosmic inert singlino background slightly colder than the cosmic neutrino background [21] • 
The phenomenology of this scenario as regards the gluino cascade decay is essentially identical to MSSM one and we 
do not make of point of trying to distinguish this type of scenario from the MSSM in this study. 

Finally we shall consider a scenario where the two lightest (predominantly inert singlino) neutralinos are both 
very light, with one around a GeV, and one much lighter, in principle in the keV range. In this case, interesting 
phcnomenoly can emerge in the gluino cascade decays as in the usual case where the lightest neutralino states are 
around half the Z mass. However, unlike that case, the lightest neutralino is not subject to direct detection limits. 
Moreover, it is possible to arrange for the correct relic abundance in such a scenario, where the lightest neutralino 
in the keV mass range is stable and constitutes Warm Dark Matter (WDM) [25]. The idea is that both the light 
neutralinos are thermally produced in the early Universe due to their couplings to the Z and Z' gauge bosons, but the 
GeV state decays late (due to its weak couplings) after both of the neutralinos have gone out of thermal equilibrium, 
and reheats the Early Universe, effectively diluting the number density of the stable keV neutralinos, such that they 
are responsible for the observed relic abundance. It is very interesting to compare the gluino cascade decays in this case 
to that where the lightest neutralinos are around half the Z mass, in order to provide an experimental "confirmation" 
of keV dark matter at the LHC. 



E. Parameter space under study 

In this study we consider a pattern of low energy soft gaugino masses that is consistent with Eq grand unification 
[7]. This typically implies that at the EWSB scale Mi « 2M± and if Mi = 150 GeV the physical gluino mass is 
around 800 GeV. In order to have a direct comparison, M\ is made equal in both the MSSM and EgSSM, 150 GeV in 
the following analysis. We also fix the physical gluino mass to be equal, 800 GeV in the following analysis. For large 
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FIG. 1. The scanned regions of the parameter spaces projected onto the plane spanned by the spin-independent cross-section, 
as i, and the relic density of the LSP, 0,-^oh 2 . The area right of the vertical solid line is excluded by WMAP |26j and 
the area above the diagonal line is excluded by XENON100, where the LSP direct detection cross-section exclusion gets 
weighted by its relic density. The 90% confidence level limit on the spin-independent LSP-nucleon cross section for a weakly 
interacting LSP with mass around 50 GeV and which makes up all of the observed amount of DM has been pushed down 
from 0.7 x 10~ 44 cm 2 = 0.7 x 10" 8 pb in 2011 [23] to 0.2 x 10" 44 cm 2 = 0.2 x 10" 8 pb in 2012 [27]. The LEP constraints on 
chargino masses (m.± > 103 GeV) and invisble Z width (F(Z -> x°X° 

) < 1.5 MeV) have been applied. The constraint on the 
Higgs mass is also taken from LEP since it holds for invisible Higgs decays, a common feature among the EsSSM points. The 
colours/shapes represent the effective gluino decay chain length ? e ff = ^'^"(O f° r each point, where P(l) is the probability for 
a chain length of I, as defined in Fig. [2] The benchmarks entitled MSSM and EgSSM-I, which are consistent with the particle 
recently discovered at the LHC being the SM-like lightest Higgs boson, are encircled. 

/i in the MSSM, and given that the effective \i in the EgSSM is large due to the limit on s coming from the limit on 
the Z' mass, there will be a neutralino that is almost the bino with a mass very close to Mi . For lower values of fx in 
the MSSM bino-Higgsino mixing occurs. 

We also consider large squark masses, with all squarks heavier than the gluino. This is motivated by the GUT con- 
strained EgSSM, where such large squark masses are a feature, although it should be noted that although the EZSSM 
scenario is consistent with such GUT constraints it has not been shown that a GUT constrained standard EgSSM 
scenario consistent with dark matter observations exists. Nonetheless we assume heavy squarks for all scenarios. 

In Fig. [l] we present the results of scans over the MSSM and the EgSSM parameter space in the (il x h 2 , a si) plane. 
The parameter space scanned over is shown in Tabs . [T] and |H| and points are linearly distributed over these ranges. The 
ranges chosen are motivated from the discussions in the previous two subsections. For both the MSSM and EgSSM 
the points are shown as long as they are consistent with the LEP limit on the SM-like Higgs boson mass, applicable 
even if the Higgs has large invisible branching fractions, but we also highlight benchmark points that are consistent 
with the particle recently discovered at the LHC being the SM-like lightest Higgs boson. For the EgSSM, where the 
LSP mass may be less than half of the Z mass, points are only shown if they are consistent with LEP limits on the 
invisible Z width, contributing less than 1-sigma. For both scans squark soft masses are set to 2 TeV, although the 
effects of squark mixing on the mass eigenstates are included. The physical gluino mass is set to 800 GeV and the 
U(1)y gaugino mass is set to 150 GeV in both cases. We do not consider scenarios in which any squarks are less 
massive than the gluino. 

We define the length of a gluino decay chain to be the number of decays after the virtual squark as in Fig. [2] We 
then also define an effective chain length for each point in parameter space 

U = X> p (0. (H.12) 

i 

where P(l) is the probability of having a decay chain of length I for that point. Intervals of effective chain length are 
colour/shape coded in Fig. [lj These scans indicate that in the EgSSM these decay chains are typically longer due to 
the bino decaying into the lower mass inert states. The distribution of effective chain lengths for these scans are also 
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TABLE I. The MSSM scanning region. A common squark 
and slepton mass scale was fixed to Ms = 2 TeV. The 
gaugino masses were fixed to Mi = 150 GeV, M2 = 
285 GeV, and M3 = 619 GeV, providing a gluino mass 
close to 800 GeV. 
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TABLE II. The E 6 SSM scanning region. A common 
squark and slepton mass scale was fixed to Mg = 2 TeV. 
The gaugino masses were fixed to Mi = 150 GeV, M[ = 
150 GeV, Mi = 300 GeV, and M- g = 800 GeV. 



plotted in Fig. 3(a) 




Decay chain length: 



FIG. 2. In the first step in a gluino decay chain the gluino decays into a quark and a squark (in our scenario it will be a 
virtual squark) which in turn decays into a second quark and a neutralino or chargino. This is the shortest possible gluino 
decay chain, which we define as having length 1 = 1. The neutralino or chargino can then decay into lighter neutralinos or 
charginos by radiating W, Z, or Higgs bosons, which typically decay into pairs of fermions. For each such decay the decay 
chain length is taken to increase by one. The radiated bosons could be on-shell or off-shell depending on the mass spectrum of 
the model. Light squarks or leptons could appear further along in the decay chain, leading to radiation of SM fermions without 
intermediate W, Z, or Higgs bosons, but in our study squarks and sleptons are heavy so this is not relevant. 



In the MSSM, for the parameters chosen, the typical decay length is 2, with the lightest chargino being initially 
produced from the virtual squark decay before subsequently itself decaying to the bino LSP. When the magnitude of 
(i is small the mixing between the gauginos (bino and wino) and the Higgsinos increases and the entire neutralino 
and chargino spectrum is pushed down. Specifically the heavier neutralinos and chargino are brought down below the 
gluino mass causing extra steps in the gluino decay chain. 

In the E 6 SSM the effective decay length is typically either 3 or 4. Initially an either charged or neutral wino is 
produced and this subsequently decays to the bino. The bino then decays into either of the two light inert neutralinos 
that are the LSP and NLSP. Which of these the bino preferentially decays into depends on the values of the 
violating couplings in blocks B a in the neutralino mass matrix in Eq. (II.7). Therefore in the E 6 SSM we typically 
expect the gluino cascade decays to be either one or two steps longer than in the MSSM. 
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left — ' ' ^"(Oi where P(l) is the probability of having a averaged over all points in the parameter scan, satisfying dark 



decay chain of length I for a point in parameter space. 



matter and collider constraints. 



FIG. 3. Statistical properties of the gluino decay chain length in the scanned parameter space. Figure |3(a)| shows how the 
effective gluino decay chain length evaluated at each point is distributed over the scanned parameter space. Most points have 
an effective decay chain length close to an integer indicating that there usually is a largely dominant decay chain length. The 
two peaks just below 3 and 4 clearly show how the EeSSM generally introduces one or two extra steps in the chain. Figure [3(b) 
show the average probabilities of different gluino decay chain lengths for both models' parameter spaces. Again, the EeSSM is 
shown to shift the probabilities to longer decay chain lengths. 



F. Benchmarks 



Below follow descriptions of the main features of the chosen benchmarks. The details of their spectrum and 
parameter values are given in Tab. Ill The diagrams for the main decay channels of the gluinos are shown in Fig. [4] 



for two main benchmarks, MSSM and EgSSM-I discussed below. The branching ratios for production of particles are 
denoted in brackets. The decay chains for both benchmarks are essentially the same up to the first two steps, with 
just a slight difference in branching ratios. The essential difference is that, in the case of the EgSSM, the lightest 
MSSM-like neutralino is no longer stable and decays in two steps to the lightest EgSSM neutralino. About 20% of 
the gluino decays go directly into the lightest MSSM-like neutralino implying a chain length I = 1 for the MSSM and 
I = 3 for the E 6 SSM. On the other hand, about 80% of the gluino decays are into a heavier neutralino or chargino, 
which subsequetly decays into the lightest MSSM-like neutralino state giving the MSSM a chain length I — 2 and 
the EgSSM I = 4. More complete diagrams showing how the gluinos decay for the different benchmarks (below) are 
shown in Fig. [15] in Appendix [A] 

The following benchmarks provide the main focus of our study: 

• MSSM: 

In this benchmark we have an LSP with a mass of 150 GeV. The low dark matter relic density is achieved via 
LSP resonance annihilation through the CP-odd and heavy Higgs bosons, A and H. The lightest Higgs boson, 
h, has a mass of 124.4 GeV and can be produced in gluino decay chains. This happens when the gluino decays 
via the next to lightest neutralino %2- The gluino decay chain length is dominantly I = 2 for this benchmark as 
for all the MSSM points scanned over. 

• E 6 SSM-I: 

In this benchmark the LSP and NLSP annihilate efficiently through the Higgs boson resonance leaving a relic 
density less than the observed relic density of dark matter. The lightest Higgs mass is around 125 GeV and the 
two lightest inert neutralino state masses are slightly above half of the Higgs mass. In this case the Higgs is 
SM-like in both its composition and its decays, since only decays into SM final states are kinematically allowed. 
If the LSP mass was be slightly below half of the Higgs mass then the Higgs boson would decay invisibly. In 
order for the two lightest inert neutralino states to be heavy enough some of the f ua f3 and fda/3 are required to be 
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MSSM: 




9 q [2i%]Xmi >Xmi Xe2 Xei 

FIG. 4. Feynman diagrams for the leading gluino decay chains for our two main benchmarks, MSSM and EeSSM-I. The 
branching ratios for production of particles are denoted in brackets. The decay chains for both benchmarks are essentially the 
same up to the first two steps, with just a slight difference in branching ratios. The essential difference is that, in the case of 
the EgSSM, the lightest MSSM-like neutralino is no longer stable and decays in two steps to the lightest EeSSM neutralino. 
About 20% of the gluino decays go directly into the lightest MSSM-like neutralino implying a chain length I = 1 for the MSSM 
and I = 3 for the EsSSM. On the other hand, about 80% of the gluino decays into a heavier neutralino or chargino, which 
subsequetly decays into the lightest MSSM-like neutralino state, giving the MSSM a chain length I — 2 and the EeSSM I = 4. 

large enough such that Yukawa coupling running becomes non-perturbative on the way up to the GUT scale. We 
estimate that here the Yukawa couplings remain purtubative up to an energy scale of order 10 12 GeV. Compared 
to the MSSM benchmark there are typically two extra steps in the decay chain as the bino-like neutralino decays 
into first the NLSP which subsequently decays into the LSP. The two extra steps in the chain make the most 
common total gluino decay chain length I — 4. The decay of the bino into the NLSP is preferred over the decay 
directly into the LSP because of the structure of the Z^ breaking trilinear Higgs Yukawa couplings given in 
Tab. ED 

• E 6 SSM-II: 

This benchmark represents a typical scenario, with decay chain features similar to in the EgSSM-I benchmark, 
where the most common decay length is 4. Yukawa coupling running remains perturbative up to the GUT scale, 
as it does in all of the following benchmarks. Just as in previous benchmark the gluino decay chain length here 
is typically two steps longer than in the MSSM because of the two extra light neutralino states, but here their 
masses are smaller. The LSP mass is below half of the Z boson mass, but decays of the Z boson into LSP-LSP 
contribute to the effective number of neutrinos as measured at LEP less than 1-sigma. Because the LSP and 
NLSP have masses not very far away from half of the Z mass they are able to annihilate relatively efficiently via 
an s-channel Z boson in the early Universe and the LSP contributes much less than the observed relic density 
of dark matter. This benchmark has a rather heavy lightest Higgs with a mass around 134 GeV. This Higgs 
decays dominantly invisibly into a pair LSPs and is ruled out if the boson candidate discovered at the LHC is 
interpreted as a SM-like Higgs boson. 

• E 6 SSM-III: 

In contrast to previous benchmarks this benchmark represents the other typical EgSSM scenario where the 
bino-like neutralino decays straight to the LSP (not via the NLSP). In this case there is only one extra step 
compared to the MSSM and the most common decay length is 3. Points with this shorter decay length are 
slightly more common when one scans over the breaking Yukawa couplings. In the same way as in E 6 SSM-II 
the LSP annihilates efficiently via an s-channel Z boson, even though it is farther away from resonance, and 
contributes much less than the observed relic density of dark matter. In this benchamark the lightest Higgs 
has a mass around 116 GeV, much lighter than in E 6 SSM-II. As above, this is ruled out if the boson candidate 
discovered at the LHC is interpreted as a SM-like Higgs boson. This benchmark represent the slightly more 
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common scenario in the parameter space where the typical decay chain has one extra step after the 150 GcV 
bino-like neutralino. 

• E 6 SSM-IV (EZSSM-I): 

This benchmark represents an EZSSM dark matter scenario as described in Ref. [53]. Here the bino-like neu- 
tralino is stable and makes up all of the observed dark matter relic density. A low enough relic density is achieved 
via the bino-like neutralino upscattering into the inert Higgsino pseudo-Dirac pair \° E3 and x%a- Since the bino 
is stable the gluino decays of this benchmark are essentially identical to those of the MSSM benchmark and also 
have the same dominant decay chain length I = 2. The bino-like neutralino's direct detection cross-section is 
small. 

This benchmark looks very similar to the MSSM benchmark and cannot be distinguished from it purely by 
analysis of gluino cascade decays. However, if a heavy or CP-odd Higgs around 300 GeV was to be excluded 
then the MSSM, but not the EZSSM, benchmark would be excluded. This relies on the fact that the EZSSM 
can have a stable bino at 150 GeV without requiring resonance annihilation through heavy Higgs boson states, 
which in contrast is required in the MSSM. The Higgs boson has a mass around 125 GeV and is SM-like in its 
composition and decays. 

• E 6 SSM-V (EZSSM-II): 

This is another benchmark of the EZSSM, but here the dark matter relic density is not explained. Here the 
inert Higgsino pseudo-Dirac pair have masses around 120 GeV, well below the bino-like neutralino. These states 
co-annihilate efficiently in the early universe and contribute a dark matter relic density less than the total 
observed dark matter relic density. The effective decay chain length is about 4. The complete decay modes for 
the gluino are shown in Fig. |15(f"J| in Appendix [Aj This benchmark provides an example of heavier lightest inert 
neutralino masses (excluding the decoupled inert singlinos) without requiring non-perturbatively large Yukawa 
couplings as in the non-EZSSM benchmark EgSSM-I. Because the lightest inert-Higgsino-like neutralinos have 
masses much larger than half of the lightest Higgs mass, which is 126 GeV, the Higgs will decay very SM-like. 

• E 6 SSM-VI (approximate Z$): 

In this benchmark the inert singlino decoupling is not exact and the inert singlino-like LSP mass is not zero, but 
has been pushed down to the 100 keV scale. This point represents the scenario where the lightest neutralino in 
the keV mass range is stable and constitutes Warm Dark Matter (WDM) [35], as described earlier, although we 
do not calculate the relic density here. 

The approximate decoupling leads to a quite long-lived bino (with a width of order 10 -11 GeV) and the step in 
the decay after the bino appears as a displaced vertex at the order of 0.1 mm from the previous step. The decay 
chain length is typically 4. The complete decay modes for the gluino is shown in Fig. 15(g) in Appendix [Aj As 



remarked, the observed relic density could be achieved by pushing the LSP mass down to the keV scale. In this 
case the last steps of the gluino decay would be likely to occur outside the detector and one would be left with 
something that looks like the MSSM. The lightest Higgs has a mass of 126 GeV and decays to the LSP are very 
suppressed due to its small mass. The Higgs could decay to the GeV scale NLSP, which in turn would decay to 
the LSP outside the detector since its width is of the order 10 -20 GeV, leading to invisible Higgs decays. The 
branching ratio for h — > X2X2 is however around 6% and is not excluded by current Higgs data. 

In the following analysis, in section IV, we use MSSM and EgSSM-I as our main benchmarks. With the exception 
of E 6 SSM-IV (EZSSM-I) the results obtained for each of the E 6 SSM benchmarks are very similar. (E 6 SSM-IV on the 
other hand looks very similar to the MSSM, since the bino is stable.) We therefore mainly give just the results for 
MSSM and E 6 SSM-I, demonstrating the qualitative difference between the MSSM and E$ models. We also include 
some results for the E 6 SSM-VI to demonstrate the effects of having an even less compact spectrum and also to show 
how little our conclusions depend on the exact spectrum of the E§ model. 



III. MODEL IMPLEMENTATION 



To scan parameter spaces of models and generate Monte Carlo events the models have to be transferred from paper 
to computer. There are various ways of doing this. Many implementations of the MSSM have already been created, 
but for the E 6 SSM there are no available sources. We have chosen to use the software package LanHEP [28] to calculate 
the Feynman rules for the EgSSM. LanHEP finds the interactions and mass mixings between the particle states in the 
model and writes an output which can be read by a Feynman diagram calculator or Monte Carlo event generator. 
LanHEP can output in several formats. We have been using the format used by CalcHEP since that is the software we 
are using for calculating Feynman diagrams and generating events. 
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TABLE III. Benchmarks motivated by the parameter scans presented in Fig. [T] and Tabs. [T] and [TTJ From top to bottom the 
classes of parameters are dimensionless input parameters; dimensionful input parameters; neutralino, chargino, and lightest 
Higgs masses (in absolute values); probabilities for certain gluino decay chain lengths; and finally dark matter properties. The 



0(±) 



are MSSM-like states, the Xm are USSM-like states, being mainly mixtures of S and B' . The Xe% are states introduced 



by the inert sector of EgSSM. The scale for squark and slepton masses is Ms = 2 TeV in all benchmarks. 



In the LanHEP implementation of the model the particle content and Lagrangian is specified. We have used a slightly 
stripped down version of the EgSSM, suitable for our purposes. What is not included from the three families of 27 
representations of the E e group are the exotic coloured states]^] their superpartners, and the inert Higgs doublets and 
SM-singlets, from the two first families. However, one should note that we are including the superpartners of these 
inert Higgs and singlet states which, as described in section [II B[ extends the neutralino and chargino sectors. The 
diagonalisation of large mass matrices appearing in the model, e.g. the 12 x 12 neutralino mass matrix, is performed 
with routines available from the SLHAplus [21] package which is well integrated with LanHEP and CalcHEP. 

The Z§ violating Yukawa couplings, x u , X4, and z, connecting the inert neutralino sector with the USSM sector, 
have been included in the LanHEP model. Turning on these couplings causes the neutralino mass matrix to leave 
its block-diagonal form and acquire non-zero off-block-diagonal elements. The model is parameterised such that the 
dimensionless input parameters are the Yukawa couplings Xijk from the SH u Hd-terms in the superpotential, the ratio 
of the Higgs doublet VEVs, tan j3, and the gauge couplings. Dimensionful input parameters of the model are the third 
generation soft trilinear scalar A-couplings, the soft masses of the squarks and sleptons, and the soft gaugino masses 



3 These can be diquarks or leptoquarks depending on the model definition. 
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-4.48 x 10 -3 


A21 


3.0 x 10" 2 


-4.34 x 1CT 2 


-8.77 x 10~ 2 





-4.34 x 10 -2 


-4.34 x 10 -2 


A12 


3.21 x 10" 2 


-3.83 x 1(T 2 


6.40 x 10" 2 





-3.83 x 10~ 2 


-3.83 x 10 -2 


An 


7.14 x 10" 4 


-1.25 x 10~ 2 


-3.36 x 10 -2 


4.22 x 10~ 2 


-1.25 x 10~ 2 


-1.25 x 10~ 2 


,fd22 


1 x 10~ 3 


-9.02 x 10~ 3 


7.06 x 10" 3 








2.0 x 10' 1 


fd21 


6.844 x KT 1 


-3.48 x 10 _1 


6.10 x 10 _1 








-3.48 x 10~ 3 


}d!2 


6.5 x 10" 1 


6.92 x 10" 1 


-7.64 x 10" 1 








6.92 x 10" 3 


fdll 


1 x KT 3 


6.17 x 10~ 3 


9.26 x 10~ 3 








6.17 x 10~ 5 


fu22 


1 x 10~ 3 


—6.77 x 10 


3.93 x 10 








1.0 x 10~ x 


}u21 


6.7 x KT 1 


-7.86 x 10- 1 


-8.56 x 10" 1 








-7.86 x 10^ 3 


/«12 


6.4 x KT 1 


2.52 x 10 _1 


-2.71 x 10 _1 








2.52 x 10 -3 


full 


1 x 1(T 3 


8.59 x 10~ 3 


-2.24 x 10" 3 








8.59 x 10" 5 


X d 2 


7.14 x 10 -4 


4.04 x 10 -3 


2.35 x 10" 4 


4.04 x 10 -5 


4.04 x 10" 5 


4.04 x 10 -5 


Xdl 


7.14 x 1(T 4 


5.11 x 10~ 4 


2.96 x 10" 4 


5.11 x 10~ 6 


5.11 x 10~ 6 


5.11 x 10~ 6 


Xu2 


7.14 x 10 -4 


-2.01 x 10~ 3 


9.04 x 10" 4 


-2.01 x 10 _s 


-2.01 x 10~ 5 


-2.01 x 10 -5 


•Eul 


7.14 x 1CT 4 


1.01 x 10~ 3 


-2.21 x 10~ 3 


1.01 x 10~ 5 


1.01 x 10~ 5 


1.01 x 10~ 5 


Z2 


1 x 10 -3 


6.02 x 10" 2 


4.16 x 10" 3 








6.02 x 10 -4 


Zl 


1 x 10 -3 


2.63 x 10 -3 


1.18 x 10~ 2 








2.63 x 10~ 5 



TABLE IV. Trilinear Higgs Yukawa couplings in the EsSSM benchmarks. The couplings A^ come from the terms XijkSiHdjH^k 
in the superpotential. Here A333 = A, A 3a/ 3 = A Q/ 3, A Q 3/3 = fdap, A a/ 33 = fnag, Mia = x da , A 3ci 3 = x ua , and A q33 = z a . 



at the electro-weak scale. One should note that we use the notation of the physical gluino mass, nig (MGo) instead of 
M3. The soft trilinear coupling A\ associated with the SH u Hd-term is exchanged for the pseudo-scalar Higgs mass 
M A . 

The details and notations of the CalcHEP model E6SSM-12.02 described above can be found in Appendix [B] and 
the model files are accessible from the High Energy Physics Model DataBase (HEPMDB) 30 . 



IV. GLUINO PRODUCTION AND DECAYS IN THE MSSM AND E 6 SSM 



The most important processes for supersymmetry searches at hadron colliders are the production of gluinos and 
squarks, provided that they are not much heavier than charginos and neutralinos. In this paper we consider the 
case where all of the squarks are heavier than the gluino, which is motivated by the GUT constrained EgSSM as 
discussed in section HE which makes the pair-production of gluinos the most attractive process for EgSSM search 
and discovery. 



A. Production cross-sections 



The tree-level cross section for gluino pair-production at the LHC at 7, 8, and 14 TeV is shown as a function of 
the gluino mass in Fifr |] The CTEQ6LL [3T] PDFs are used and the cross section is evaluated at the QCD scales 
Q = Vs and Q = nig rl 

In Fig. [5] one can see a large scale dependence of the cross section due to the uncertainty in the leading order 
calculation, which is substantially reduced at NLO level [32 35 . At Q = m g , for which the cross section is about 50% 
larger than at Q = y/s at tree level, the product of NLO and NLL K-factors is in the 2.5-5 range for y/s = 7 TeV for 
the 500-1500 GeV mass range [35]. To be on the conservative side we use LO cross sections evaluated at Q = nig in 
our analysis. For our benchmarks with a gluino mass of 800 GeV, the production cross sections are 20.6 fb, 47.5 fb, 
and 839 fb for the y/s — 7,8 and 14 TeV respectively. So, at the present y/s — 8 TeV LHC energy and expected 
integrated luminosity of 20 fb -1 by the end of 2012 about 1000 gluino pairs would have been produced. 



4 Both renormalisation and factorisation scales were chosen to be equal. 
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FIG. 5. The tree-level cross section for gluino pair production as a function of the gluino mass, m§. The solid (or dashed) 
lines represent, from bottom to top, the LHC at 7 TeV (black), 8 TeV (red), and 14 TeV (blue). The CTEQ6LL set is used 
for PDFs. The QCD scale, Q, is set to the gluino mass, Q = m§, for the solid lines and to the centre of mass energy, Q = 
for the dashed lines. The scale dependence of the cross section is an effect of the uncertainty of the leading order calculation. 
Including NLO corrections is known to bring up the cross section by at least a factor 2 [32], so we are underestimating the 
production rate for gluinos slightly with this leading order calculation. 



B. Signatures and distributions 



Since the EgSSM introduces new neutralinos, naturally lighter than the MSSM-Uke LSP, the gluino decay chains 
will be longer than in the MSSM in general. This is confirmed and illustrated by the parameter scans in Fig. [I] and the 
benchmarks defined in Tab. |III| In order to study the LHC phenomenology of gluino cascade decays of the EgSSM, and 
the MSSM for comparison, we have performed Monte Carlo analyses using the CalcHEP [53] package with CTEQ6LL 



[5T] PDFs. With the exception of the multi-jet analysis in section IV Dl we restrict ourselves to a parton-level 



analysis. We do however take into account a realistic electromagnetic energy resolution, given by 0.15/ v /^(GeV), 
typical for the ATLAS and CMS detectors, as well as their typical hadronic energy resolution of 0.5/^/E{GeV) and 
perform the respective Gaussian smearing for leptons and quarks. We define leptons (jets) by requiring px > 10 GeV 
(20 GeV) and \rj\ < 2.5 (4.5) and a lepton isolation of Ai?(lepton, jet) > 0.5. 

The longer decay chains of gluinos lead to less missing momentum, p™ lss , and larger effective mass M c g — p'^ lss + 
S visible I.PT Slble l> as measured in the detector as one can see in Fig. which presents the respective distributions for 
our main benchmarks (MSSM and EgSSM-I). Although the p™ lss distribution is quite different, one should note that 
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the effective mass distribution is not significantly different between the models. This happens because the effect from 
the suppressed missing momentum in the case of the EgSSM is partially canceled by the effect of the increase of visible 
momentum, due to the longer gluino cascade decay. There is a slight overall increase of the effective mass due to the 
fact that visible momenta are added up as magnitudes while the missing momentum is a vectorial sum. The reduced 
amount of missing transverse momentum in the EgSSM makes it less discoverable, compared to the MSSM, in typical 
SUSY searches which focus on all-hadronic events with large missing momentum. 

Another important feature of the long decay chains of the EgSSM is the increase in lepton as well as jet multiplicity, 
as shown in Fig. [7J again for the benchmarks MSSM and EgSSM-I. This feature allows us to rely on multi-lepton 
requirements for background reduction rather than cuts on missing energy. There is a significant loss of statistics by 
using this strategy, however it turns out to be a very important channel for discovery of gluinos with long decay chains 
and indeed a channel in which the EgSSM is largely dominant compared to the MSSM. This makes the multi-lepton 
channels essential for distinguishing the models. 




FIG. 6. Missing transverse momentum (left) and the effective mass (right) before cuts for the MSSM and EeSSM-I benchmark 
with trig — 800 GeV at \/s = 7 TeV. The EgSSM predicts significantly less missing transverse momentum and slightly larger 
effective mass compared to the MSSM. The longer gluino decay chains of the EeSSM, with a lighter LSP in the end, provide 
less missing and more visible transverse momentum. The effective mass does not distinguish the features of these models since 
it is a sum of visible and missing transverse momenta. 



C. Searches at y/s = 7 TeV LHC 

There has not been any indications of SUSY from the LHC during its run at s/s = 7 TeV. We have investigated 
different SUSY search channels at this energy to understand the status of our benchmarks and what limits can be put 
on the EgSSM and which channels we expect to be the most favourable for discovery and distinguishing the models. 
We compare our signals with published backgrounds used by CMS and ATLAS at this energy. We have scaled all the 
channels to an integrated luminosity of 10 fb _1 for comparison, which is approximately the amount of 7 TeV data 
acquired by the two experiments. Benchmarks with an 800 GeV gluino mass are considered here. 

1. leptons 

The long gluino cascade decays with less missing momentum would be less visible in the main SUSY searches 
based on jets and missing energy (see e.g. [2] and [37]) which provide the best statistics and strongest exclusions for 
MSSM. In these searches the EgSSM parameter space is less constrained as compared to the MSSM and the acquired 
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FIG. 7. Lepton multiplicity (left) and jet multiplicity (right), requiring pr > 10 GeV, \rj\ < 2.5, and Ai?(lepton, jet) > 0.5 
for l epto ns and pr > 20 GeV and \rj\ < 4.5 for jets. The benchmarks considered are the MSSM and EeSSM-I as presented in 
with rrig — 800 GeV. The LHC setup is used with ^fs = 7 TeV and normalised to 10 fb _1 of integrated luminosity. Due 



III 



Tab. 

to the longer gluino decay chains of the EgSSM it predicts many more visible particles in collider experiments, both leptons 
and jets. This suggests that searches for EgSSM gluinos should be more favourable in multi-lepton and multi-jet searches. 



exclusions do not hold for this model. The main reason for this is the hard cuts on missing energy and its ratio to 
the effective mass since the distributions for these variables are significantly different for MSSM versus EgSSM as we 
demonstrate here. 

The effective mass distribution for our benchmarks for an 800 GeV gluino mass is plotted on the top of the 
backgrounds from ATLAS and CMS in Figs. |8(a) and |8(b)| after all cuts have been applied except for the final 
selection cut on the effective mass itself. The signal from EgSSM is suppressed as compared to the MSSM and, more 
importantly, both benchmarks are well below the background, illustrating the difficulty of discovering SUSY at the 7 
TeV LHC in the case where the gluino mass is around 800 GeV, assuming the squarks to be much heavier. 

Even though the 0-lepton signature with a jet multiplicity of about four is not favoured for EgSSM hunting, the 0- 
lepton signature for this model could be still interesting for the cases of larger jet-multiplicity. For multi-jet ch annels, 
analyses beyond the parton level are essential. We discuss this in detail for the case of y/s = 8 TeV in section IV D 1 



where we perform one example of a beyond-the-parton level analysis. Apart from this, we restrict ourselves here by 
the parton level analysis and in particular we shall focus on the tri-lepton signature. 



2. 1-2 leptons 



The selection of events with leptons provides easy triggering and efficient background suppression at the cost of 
worse statistics. To exemplify this we compare the signal distributions for the benchmarks versus the backgrounds 
for two ATLAS searches, a single lepton search in Fig. 8(c) and a two same-sign lepton search in Fig. |8(d)"| One can 



see that the signal-to-background ratio is better in these leptonic searches compared to the all-hadronic searches. In 
the effective mass distribution of the 1 muon channel from ATLAS's 1 lepton plus 4 jets search shown in Fig. |8(c " 



one 



can see that the signal level is not extremely far below the background level in the high effective mass region. The 
EgSSM signal is still suppressed as compared to the MSSM in thi s one lepton search, however, considering the p™ ISS 
distribution for the two same-sign lepton search by ATLAS in Fig. |8(d)] one sees how the EgSSM signal overtakes the 
MSSM benchmark's by requiring one more lepton. This is due to the fact that two same-sign leptons in the final state 
become more likely in the EgSSM than in the MSSM, simply because it predicts more leptons in general. Even though 
the EgSSM signal has got stronger than the MSSM signal in this 2SS channel compared to the one lepton channel, 
the signal-to-background ratio now looks worse. This is because the choice of using the missing momentum instead 
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FIG. 8. Distributions for published 0-4 lepton searches at 7 TeV, scaled to 10 fb~ for comparison. The signal contributions 
from the MSSM and EsSSM-I benchmarks with m, = 800 GeV are plotted on top of published backgrounds. The distributions 
shown are after all cuts except the final selection cut on the plotted distribution. lepton searches (with jets<4) give bad 
background suppression and favour MSSM. Requiring two or more leptons makes MSSM more suppressed. For multi-lepton 
searches the signal to background ratio for the EeSSM is better but the signal statistics is low. 



of the effective mass as the variable to define the signal region is not favourable for the EgSSM since the respective 
signal is very vulnerable to hard cuts on this variable. Also for higher multiplicity searches, using the effective mass 
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to define the signal region allows for the signal-to-background ratio to be improved as we will show below. 



3. 3-4 leptons 



Requiring additional leptons makes the statistics even worse, but it allows the signal to appear above the background 
enough to allow a reasonable signal significance in order to test the models under study. Comparison of the benchmark 
signals with the background from a 3 lepton search by CMS is shown in Fig. 
in Fig. |8(f)| The three lepton channel defined by the cuts 



(e) and from a 4 lepton search by ATLAS 



pr(h) > 20 GeV 

p T (h) > 10 GeV (IV.l) 
p T (h) > 10 GeV 



where I = fi or e with < 2.5, and Ai?(lepton, jet) > 0.5, seems promising, showing a possible excess in the 
high missing transverse momentum region. The background used by CMS is not evaluated for large enough missing 
transverse momentum however. To explore the large missing transverse momentum region and to expand the analysis 
further we have produced backgrounds for this channel, including several processes, using CalcHEP. 

The dominant backgrounds come from ZWj and ttV . Other important contributions come from ZW and tt. 
We also considered backgrounds such as ZWjj, ZZ and ZZZ which we found to be subleading. Our background 
predictions at y/a = 7 TeV agree well with backgrounds used in the multi- lepton searches by CMS [40] and ATLAS [41] . 
They only major difference is in the very low end of the p™ lss distribution where the CalcHEP generated backgrounds 
are suppressed. 

This difference in the transverse missing momentum distribution between our results and the results from the full 
detector simulation occurs because we do not simulate any source of instrumental missing energy in our analysis. 
However, this difference does not affect the our results since we are not using missing pt information directly, similar 
to approach of [35], and, moreover, this difference effectively vanishes after the cut is applied on the effective mass 
variable as used at last stage of our analysis. 

We would like to stress that parton level analysis are quite accurate to the tri-lepton signature we study in this 
paper. Signatures with lower lepton multiplicity require considering QCD backgrounds from jets faking leptons and 
instrumental defects which could affect the low p™ lss region. This background is difficult to take into account at the 
parton-level event generator and even at the level fast detector simulation. Therefore analysis of signatures with the 
lepton multiplicites below three are outside of the scope of our paper. In case we compare signal versus background 
for these signatures for the illustration purposes we chose to rely on published backgrounds for those channels in this 
paper. 



The result for the p™ lss distribution is shown in Fig. 9(a) One can see that the EgSSM signal is now at the same 



level as the background and maybe a little higher for large p™ lss . If one instead considers the effective mass as a 
selection variable for the three lepton channel, the situation looks much more promising, at least for the EgSSM. This 
can be seen from the effective mass distribution presented in Fig. [9(b)] CMS has not been using the effective mass 
to define the signal region for this channel but uses the missing transverse energy and the hadronic transverse energy 
instead. Our way of defining the signal region by the effective mass is on the other hand much closer to the way 
presented by ATLAS in J3TJ or as suggested in [42]. A cut on M eff at 950 GeV gives S = 11.5 signal events for the 800 
GeV gluino mass E 6 SSM-I benchmark and B — 0.4 background events, providing an expected 5.7a excess at 10 fb _1 , 
using the definition of statistical significance S12 = 2(\/5 + B — s/B) valid for small statistics [351 133] . 
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FIG. 9. Distributions for p™ lss (a) and M off (b) at the LHC with 10 fb^ 1 at y^i = 7 TeV after requiring at least three leptons. 
The benchmarks both have a gluino mass rrig — 800 GeV. Backgrounds have been generated by CalcHEP and agree well with 
published ones such the ones by CMS, shown in Fig. |8(e)] The EeSSM-I benchmark is shown to present a signal larger than the 
background for large missing momentum (a) even though it is a model that predicts quite small amounts of missing momentum. 
The signal presents itself more strongly in the effective mass variable where there is no need for a cut on the missing transverse 
momentum. 



D. Searches at s/s = 8 TeV LHC 

1 . 6 jets 

Since the multi-jet channels could provide good prospects of discovery of and differentiation between benchmarks 
we investigate the effect of the cuts, 



rrnmss 
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Pr(h) 
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> 0.25(6j) 


m eS 


> 1300 GeV, 



(IV.2) 



used by ATLAS [45) for 6-jets analysis applied to our benchmarks. To perform this kind of multi-jet analysis we 
are forced to go beyond our parton-level analysis. The need for a more dedicated analysis, with hadronisation and 
detector effects, comes mainly from the importance of initial and final state radiation that plays an important role 
when requiring more than four jets in the gluino decays. This is because the first four jets are well approximated by 
our parton level analysis since they typically are the ones that originate straight from the gluino decay. To generate 
the events at the level of fast detector simulation we fed events from CalcHEP to the PGS 46J package. We also 
compare the signal from the MSSM and E 6 SSM with an mSUGRA point which is excluded, but very close to the 
limit. The result for the effective mass distribution for the signals before and after cuts are shown in Fig. [10] 
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Our result, produced with CalcHEP and PGS, is in good agreement with the ATLAS result and is consistent with 
the experimental data. 
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FIG. 10. Comparison between mSUGRA, MSSM, and EeSSM benchmarks in the 6 jet channel, E-medium, used by ATLAS 
[45] . The effective mass distribution for the gluino signal are plotted before (a) and after (b) cuts at ^/s — 8 LeV and 10 fb~ 
of integrated luminosity. The events left after the last signal region cut on the effective mass, M e g > 1300 GeV, are given in 
Tab. [V] The benchmarks with solid lines all have a gluino mass of 1 TeV while the benchmarks with dashed lines have a gluino 
mass of 800 GeV. After cuts (b) the EgSSM benchmark with an 800 GeV gluino mass has a distribution not very different from 
the mSUGRA point with a 1 TeV gluino mass. 



Our analysis shows that, also in the multi-jet channels, the EgSSM will be suppressed compared to the MSSM 
and mSUGRA models, mostly because of its small missing transverse momentum. The larger jet multiplicity of the 
EgSSM at the parton level turns to be not such a great discriminator between two models, at least for this signature, 
because the initial and final state radiation allows most of the signal from the MSSM model to pass the 6-jets selection 
requirement. Even though it might still be a good discovery channel for different selection/analysis for the EgSSM, 
we note that, for the current selection, the EgSSM signal from an 800 GeV gluino is about at the same level as the 
signal from the mSUGRA model with a 1 TeV gluino. This is the very important point we would like to convey in 
this paper: the MSSM limits are not quite applicable to the EgSSM, for example the gluino mass limit could easily 
differ by about 200 GeV between two models in the heavy squark limit as we demonstrate above. Therefore the status 
of the EgSSM is quite different even from generic MSSM one for the current LHC analysis, which needs to be tuned 
in order to explore EgSSM parameter space. 



mj/TeV Events 



mSUGRA (1850,380) 


1 


6.18 


MSSM 


1 

0.8 


5.59 
18.16 


E 6 SSM 


1 

0.8 


1.05 
5.58 



TABLE V. The events left at 8 TeV and 10 fb _1 for benchmarks from three models after the E-medium set of cuts, including 
the final cut on the effective mass, M e g > 1300 GeV. Here the EgSSM benchmark with a 800 GeV gluino mass is left with less 
events than the mSUGRA benchmark with a 1 TeV gluino mass. 
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2. 3 leptons 



Performing the same analysis for the promising 3 lepton channel at the LHC with ^fs — 8 TeV as was done for 
the case with y/s — 7 TeV we are able to calculate the discovery prospects for our benchmarks. Again, the 3 lepton 
signature is defined by the cuts in |IV.l| In Fig. [IT] the effective mass distributions are plotted for our benchmarks and 
SM backgrounds for 20 fb _1 of integrated luminosity. The different plots show how the effective mass distributions 
change for the benchmarks if the gluino mass is varied. 
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FIG. 11. Plots of M ett after requiring at least 3 leptons at LHC at *Js =8 TeV. The integrated luminosity is taken to be 20 fb _1 . 
The different subfigures show the signal distributions for the MSSM, EgSSM-I, and EeSSM-VI benchmarks for different values 
of the gluino mass. The EgSSM-VI benchmark is similar to E6SSM-I, but the lighter LSP mass, and thus larger mass gap 
between it and the bino-like neutralino, causes the signal to be stronger since higher pr leptons are more likely to be produced. 
The distributions for gluino masses of 700 GeV (a), 800 GeV (b), 900 GeV (c), and 1000 GeV (d) show that the signal to 
background ratio is not affected much as the gluino mass increases, but the statistics become bad since the cross section gets 
small. 



In these plots we have included another benchmark, E 6 SSM-VI, for comparison with the E 6 SSM-I and MSSM 
benchmarks. This benchmark gives a signal slightly above that of the EgSSM-I. This is because the EgSSM-VI 
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benchmark has a long gluino decay chain with an even less compact spectrum. The larger mass difference between 
the MSSM-like lightest neutralino and the inert singlinos implies that higher pt leptons can be radiated from that 
step in the decay chain. This causes an increase in the number of lepton surviving the lepton identification cuts. 

A final cut, defining the signal region, is made on the effective mass. We let the signal region cut depend on the 
gluino mass to enhance the expected significance and define it as M e R > lArrig. 



Using the definition of statistical significance, S%2 — 2(V S + B — \AB), valid for small statistics we calculate 

the expected excess for different gluino masses using our mass dependent signal region cut. 

The significance is plotted as a function of the gluino mass in Fig. 1 2 (a) | where a K- factor of 3 has been applied to 
the signal. The expected number of events for the EgSSM-I and EgSSM-VI benchmarks with gluino mass nig = 900 
GeV and the background before and after the final cut on the effective mass (M c g > \Am~ g = 1260 GeV) are presented 
in Tab. |VI| The table also lists the expected significances with and without a K-factor of 3 for the two benchmarks. 
The integrated luminosity needed for discovery and exclusion of a particular gluino mass in the E 6 SSM in the 3 lepton 



channel is shown in Fig. 12 where again a K-factor of 3 has been applied to the signal. The plot shows that a 2a 
exclusion of gluino masses below 1 TeV is possible with data acquired by ATLAS or CMS at the end of the year 2012. 
The MSSM benchmark is still well below the background at this stage and is therefore not included in these plots. 
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TABLE VI. The expected number of events after the first and second cuts in the 3 lepton analysis for the EeSSM-I and 
EsSSM-VI benchmarks with rrig = 900 GeV and SM background at 20 fb _1 and 8 TeV. Also the significance, based on signal 
and background events after the second cut, with K-factors of 1 and 3 applied to the signal. 
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FIG. 12. The gluino mass reach at yfs =8 TeV for the three lepton channel. The gluino mass is varied for the benchmarks 
EsSSM-I and EeSSM-VI to estimate the expected significance for different gluino masses. The significance is calculated with the 
events remaining after a selection cut requiring M c g > 1.4m§. A K-factor of 3 has been applied to the signal. The EgSSM-VI 
benchmark (shown in blue) is more accessible for exclusion or discovery than benchmark EeSSM-I since it has a bigger mass 
gap between the bino-like and inert-singlino-like neutralinos, providing higher pr leptons. 



E. 3 lepton searches at yfs = 14 TeV LHC 



At higher collider energy the cross section for gluino production increases considerably. This causes both our MSSM 
and EgSSM benchmarks to be clearly visible above the background, as can be seen in Fig. [13] The figure shows the 
effective mass distribution for the MSSM, E 6 SSM-I, and E 6 SSM-VI benchmarks for different gluino masses in different 
subfigures, where a requirement of at least three leptons has been applied in the same way as in the 8 TeV analysis. 
For all gluino masses the EgSSM-VI benchmark gives the largest signal, just as in the 8 TeV scenario discussed in 
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The 14 TeV collider at such a large integrated luminosity as 100 fb , which is used in Fig. 13 allows 
for statistics needed for discovery of high mass gluinos. Heavier gluinos would push the effective mass distribution to 
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higher values where there is essentially no background and where it is just a matter of acquiring enough statistics, 
but once that is done one expects a very clean signal. 
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FIG. 13. These plots show the evolution of the effective mass distribution after requiring at least three leptons for the 
benchmarks MSSM, EeSSM-I, and EeSSM-VI as the gluino mass changes. In Figure 13(a) the gluino mass is 900 GeV and 
at this large integrated luminosity of 100 fb _1 at */s = 14 TeV the MSSM benchmark is almost discoverable. The EgSSM 
benchmarks are both well above the background and clear signals are expected due to the good statistics. 



In Fig.[l4]the gluino reach for the benchmarks at the 14 TeV collider is presented in the same way as was done for the 
8 TeV collider in Fig. [12] again with a K-factor of 3 applied to the signal. The only difference is that the gluino mass 
dependent cut on the effective mass which defines the signal region is taken to be M e g > nig instead of M e s > lArrig. 
For the MSSM benchmark, which has become accessible at this energy, one will be able to exclude gluino masses up 
to ^1400 GeV through this channel for an integrated luminosity of 100 fb _1 . The expec ted sig nificance for the MSSM 
is however about an order of magnitude below the EgSSM benchmarks as shown in Fig. |14(a) An order of magnitude 
difference in significance implies two orders of magnitude differen ce for the the integrated luminocity required for 
exclusion of a particular gluino mass, as can be seen in Fig. 14(b) For the EgSSM-I and EgSSM- VI the 100 fb _1 of 
data at 14 TeV allows not only for exclusion, but potentially for a 5<r discovery for the whole gluino mass range up 
to 1.5 TeV. 
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FIG. 14. The gluino mass reach at ^J~s =14 TeV for the three lepton channel. The gluino mass has been varied for the 
benchmarks EgSSM-I, EeSSM-VI, and MSSM. The expected significance is calculated using the signal and background events 
surviving the cut M c h > rrig. A K-factor of 3 has also been applied to the signal. Here at the higher collider energy the MSSM 
benchmark also becomes discoverable through this 3 lepton channel. The EgSSM benchmarks are discoverable up to almost a 
1500 GeV gluino mass. 

It should be stressed once more that the results for the MSSM benchmark are different from other studies since we 
are not dealing with a GUT constrained model, but an electro-weak scale model and the spectrum is more general. 
Our study focuses on a specific scenario in which the squarks are two to three times heavier than the gluino. Therefore, 
typical results obtained in previous multi-lepton analyses for SUSY searches differ from our results. As an example, the 
three lepton signals derived for the mSUGRA points in [32] are larger than the signals for the benchmarks considered 
in this paper with the same gluino mass. This is mainly due to the lighter squark masses assumed in that study, 
giving rise to a larger cross section, but also partly due to the entire assumed spectrum being very different. These 
differences plays a crucial role in providing different pt distributions for the leptons arising from gluino cascade decays 
between the different models. 



SUSY is an attractive candidate for new BSM physics, but so far it remains elusive at the LHC. A generic prediction 
of all SUSY models is the gluino, which can be produced with a large cross-section at the LHC due to its colour octet 
nature. Moreover in many SUSY models it may be the lightest coloured SUSY particle, which could make it the first 
SUSY particle to be discovered. However the gluino typically decays in cascade decays, via a chain of charginos and 
neutralinos, emitting jets or leptons at each stage of the decay chain, with the chain ending when the LSP is produced 
which is typically the lightest neutralino. 

In this paper we have discussed the gluino cascade decays in Eq models, which include matter and Higgs filling 
out three complete 27 dimensional representations close to the TeV scale. In such models there are three families of 
Higgs doublets of both H u and Hd kinds, plus three Higgs SM-singlets of the NMSSM kind S. Only the third family 
acquire VEVs and the other two families do not and are called inert. All these Higgs states are accompanied by 
spin-1/2 SUSY partners, the Higgsinos. The extra inert charged and neutral Higgsinos, some of which are necessarily 
light, will mix with the usual MSSM charginos and neutralinos, yielding extra light states, providing extra links in 
the gluino cascade decays, and hence extra jets and leptons, with less missing energy. 

The extra neutralinos and charginos generically appearing in a large class of Eq inspired models lead to distinctive 
signatures from gluino cascade decays in comparison to those from the MSSM. These signatures involve longer decay 
chains, more visible transverse energy, higher multiplicities of jets and leptons, and less missing transverse energy 
than in the MSSM. We have studied this effect in gluino cascade decays for the MSSM and E^SSM and have shown 
that it can provide a characteristic and distinctive signature of the model, enabling an early discovery of the -E^SSM 
which may be discriminated from the MSSM. 

In order to demonstrate this, we have first defined a rather large set of benchmark points within the EgSSM. These 
benchmark points are chosen for the variety of ways in which the inert charginos and neutralinos can appear, and 
they all are chosen to have a nominal gluino mass of about 800 GeV, although this can be varied while keeping the 
inert chargino and neutralino masses unchanged. The results discussed below are remarkably robust, and apply to all 
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of the Eq benchmark points. These Eq benchmarks are also compared to an MSSM benchmark which is chosen to 
have similar conventional (non-inert) chargino and neutralino masses to the Eq ones, as well as mSUGRA points, in 
order to verify the model independence of the conclusions. 

Given this set of benchmark points, we have then studied gluino pair production and decay at the LHC, first at 
7 TeV, then at 8 TeV, and eventually at 14 TeV using a Monte Carlo analysis. We already know that the gluino was 
not discovered at 7 TeV, which motivates the nominal choice of gluino mass of 800 GeV. For this gluino mass, we have 
first studied the missing transverse momentum and effective mass distributions for representative benchmark points, 
and seen that the former has a softer spectrum in the Eq models when compared to the MSSM, as expected, while the 
latter has a similar or slightly harder distribution in the effective mass variable. Staying at 7 TeV, we then calculated 
the lepton and jet multiplicities in Eq models and showed that they are significantly higher than in the MSSM. This 
motivated a study of lepton channels with increasing numbers of leptons, and decreasing statistics, where we showed 
that the 3 lepton channel provided a very good discriminator between the MSSM and the Eq models, although the 
statistics are rather too low at 7 TeV. 

At 8 TeV we studied the 6 jet channel and showed that, in an effective mass analysis, a 1 TeV MSSM gluino may 
provide a similar limit to that of an 800 GeV Eq gluino. Turning to the promising 3 lepton channels at 8 TeV, we find 
increased statistics and possible observable signals in this channel for a range of gluino masses, which would provide 
a striking confirmation of Eq models. We calculate the required integrated luminosity in order to either discover or 
exclude the Eq gluino in this channel at 8 TeV. We emphasise again that the MSSM gluino is unobservable in the 
3 lepton channel. Finally we have repeated the analysis for the promising 3 lepton channel at 14 TeV and found 
analogous results for required integrated luminosity to discover or exclude the Eq gluino there as well. 

In conclusion, the Eq inspired models are clearly distinguishable from the MSSM in gluino cascade decays at 
the LHC with the full data set at 8 TeV, and certainly at 14 TeV, using the three lepton channel that we have 
proposed. Moreover, the present limits on the gluino mass, for example from a multijet analysis, are weaker (and 
therefore not applicable) for Eq models in comparison with the MSSM, due to the longer decay chains with less 
missing transverse energy that is expected in Eq models. We emphasise to the LHC experimental groups that the 
distinctive features present in gluino cascade decays, resulting in different search strategies including the choice of the 
kinematical variables, such as those discussed here, not only represents a unique footprint of a particular model but 
also may provide the key to an earlier discovery of supersymmetry. 
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Appendix A: Decay Diagrams for Benchmarks 



This section displays the possible gluino decays of each benchmark in a more complete fashion. Horizontal lines 
representing particle states are separated proportionally to their mass difference. Some exceptions are made where 
states are closely degenerate, in which case the lines have been separated more. An arrow then connects lines 
where possible decays occur with the corresponding branching ratio written above. The gluino decays of the MSSM 
bench mark is shown in Fig. [15(a)] E 6 S SM-I in Fig. |15(b)[ E 6 SSM-I I in Fig. [15(c)) E 6 SSM-III in Fig. [l5(d)| E 6 SSM-IV 



in Fig. 15(e) E 6 SSM-V in Fig. 15(f) and E 6 SSM-VI in Fig. 15(g] 
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FIG. 15. Gluino decay diagrams for the MSSM and EeSSM benchmarks, showing the leading decays (contributing more than 
90%) for the involved sparticles. The vertical line spacings are proportional to the mass splitting among the particles. 
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Appendix B: Details about the CalcHEP Model E6SSM-12.02 



In this section the contents and properties of the CalcHEP model E6SSM-12 . 02 (hcpmdb:1112.0106)[T7] are described 
in detail. The model files are accesible from HEPMDB [3U] where one can either run it with CalcHEP on the IRIDIS 
cluster via the web interface or download it and run it on one's local CalcHEP installation. The model constitutes 
of four files; varsNN.mdl, prtclsNN .mdl, funcNN.mdl, and lgrngNN.mdl; for input variables, particle definitions, 
functions or constraints, and Feynman rules for vertices, respectively. 



1. Particle content 



The model shares many features and particles with the MSSM and its extensions, e.g. the USSM or NMSSM, but 
has a greater particle content. This CalcHEP version of the EgSSM, E6SSM-12.02, includes particles from the three 
generations of 27 representations but not all. What is yet to be implemented is the full Higgs sector, including the 
inert Higgs boson states, and the coloured exotics. This is work in progress and will be added in a later version of the 
model. The extra particles included compared to the MSSM are a SM-singlet S, which mixes with the two MSSM-like 
CP-even Higgs particles; a Z'-boson; two extra chargino states from inert Higgsinos; and 8 extra neutralino states, two 
from the bino' and the singlino and six from the inert neutral Higgsinos and singlinos. The CalcHEP particle names 
and properties as used in the prtclsNN .mdl file are shown in Tab. VII The table shows the full name of each particle 
followed by the symbols used in CalcHEP for it and its antiparticle. The PDG code is a positive number assigned to 
each particle for referencing in the code, the antiparticles have PDG codes equal to the negative of the particles' PDGs. 
The spin, mass, width and color properties are also given in CalcHEP notation. Following the convention from earlier 
SUSY models in CalcHEP, superparters, i.e. particles which are odd under i?-parity, are denoted with a - prefix. The 
conventions we use for the particle mass and width variables are M< particle symbol > and w< particle symbol >. If 
the particle is charged the particle symbol is modified with a suffix -c for "charged" to avoid confusion, e.g. MHc for 
the charged Higgs H+. In the case of superpartners the ~ is not used in the mass or width variable. Instead a prefix 
S- or suffix -o is attached to the particle name depending on whether the sparticle is a Sfcrmion or a bosono, e.g. 
MGo for the gluino mass and MSel for the first selectron mass. A number is added to the variable name whenever 
it refers to particles for which there are more than one of the same type, e.g. neutralinos or stops. We try to keep 
as close as possible to these conventions without departing too much from conventions used in earlier SUSY models. 
There are still some deviations from the set up conventions, e.g. the charged slepton sector, which we hope will not 
be too confusing until a more uniform way of naming the particles and variables is introduced. In future versions of 
the model the conventions will be used more strictly. 



2. Input parameters 



The input parameters used in the model, with the exception of some SM parameters, are listed in Tab. VIII The 
pseudoscalar Higgs mass is used as an input variable instead of the soft trilinear lambda coupling from the term 
A\\SH<iH u . The electroweak soft gaugino masses Mi, M2, and M[ are denoted MG1, MG2, and MGlb, respectively. 
The physical gluino mass is denoted MGo. The soft squark masses for the first two generations are set by a common 
squark mass scale, MSq, by default. It is however easy to modify the model files by moving the soft squark masses 
Mql, Mq2, Mul , Mu2, Mdl, Md2 from the function file funcNN.mdl to the input parameter file varsNN.mdl and 
assigning them separate values. 



3. Functions and dependent parameters 



In the file funcNN.mdl all dependent parameters are listed in terms of input parameters and dependent parameters 
above themselves in the list. The dependence can be given as simple algebraic expressions or as functions of external 
functions. As an example, expressions for mass matrix elements are given in this file. These matrix elements are 
then used as inputs for numerical diagonalisation routines from the SLHAplus library, which comes with CalcHEP. 
SLHAplus then returns evaluated particle masses, which arc defined in this file. As an example of what the contents 
of this file look likes, a few lines from it are presented in Tab. |IX| 
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Full Name 


particle 


antiparticle PDG ID 2 x spin 


mass 


width 


c 


gluon 


G 


G 


21 


2 








8 


neutrino 


nl 


Nl 


12 


1 








1 


electron 


el 


El 


11 


1 


Mel 





1 


muon-ncutrino 


n2 


N2 


14 


1 








1 


muon 


e2 


E2 


13 


1 


Me2 





1 


tau-neutrino 


n3 


N3 


16 


1 








1 


tau-lepton 


e3 


E3 


15 


1 


Me3 





1 


u-quark 


u 


U 


2 


1 


Mu 





3 


d- quark 


d 


D 


1 


1 


Md 





3 


c-quark 


c 


C 


4 


1 


Mc 





3 


s-quark 


s 


S 


3 


1 


Ms 





3 


t-quark 


t 


T 


6 


1 


Mt 


!wt 


3 


b- quark 


b 


B 


5 


1 


Mb 





3 


light Higgs 


hi 


hi 


25 





Mhl 


!whl 


1 


heavier Higgs 


h2 


h2 


26 





Mh2 


!wh2 


1 


heaviest Higgs 


h3 


h3 


27 





Mh3 


!wh3 


1 


pseudoscalar Higj 


;s ha 


ha 


28 





Mha 


!wha 


1 


charged Higgs 


H+ 


H- 


37 





MHc 


!wHc 


1 


photon 


A 


A 


22 


2 








1 


Z-boson 


Z 


Z 


23 


2 


MZ 


!wZ 


1 


W-boson 


w+ 


W- 


24 


2 


MW 


!wW 


1 


Z-primed-boson 


Zb 


Zb 


32 


2 


MZb 


!wZb 


1 


chargino 1 


-1+ 


-1- 


1000024 


1 


MCol 


IwCol 


1 


chargino 2 


-2+ 


-2- 


1000037 


1 


MCo2 


!wCo2 


1 


chargino 3 


-3+ 


-3- 


1000038 


1 


MCo3 


!wCo3 


1 


chargino 4 


-4+ 


-4- 


1000039 


1 


MCo4 


!wCo4 


1 


neutralino 1 


~ol 


~ol 


1000022 


1 


MNol 





1 


neutralino 2 


~o2 


~o2 


1000023 


1 


MNo2 


!wNo2 


1 


neutralino 3 


~o3 


~o3 


1000025 


1 


MNo3 


!wNo3 


1 


neutralino 4 


~o4 


~o4 


1000026 


1 


MNo4 


!wNo4 


1 


neutralino 5 


~o5 


~o5 


1000027 


1 


MNo5 


!wNo5 


1 


neutralino 6 


~o6 


~o6 


1000028 


1 


MNo6 


!wNo6 


1 


neutralino 7 


~o7 


~o7 


1000029 


1 


MNo7 


!wNo7 


1 


neutralino 8 


~o8 


~o8 


1000030 


1 


MNo8 


!wNo8 


1 


neutralino 9 


~o9 


-o9 


1000031 


1 


MNo9 


!wNo9 


1 


neutralino 10 


~oA 


-oA 


1000032 


1 


MNoA 


IwNoA 


1 


neutralino 11 


~oB 


-oB 


1000033 


1 


MNoB 


IwNoB 


1 


neutralino 12 


~oC 


-oC 


1000034 


1 


MNoC 


IwNoC 


1 


gluino 


~g 


~g 


1000021 


1 


MGo 


!wGo 


8 


1st selectron 


-el 


-El 


1000011 





MSel 


IwSel 


1 


2nd selectron 


~e4 


-E4 


2000011 





MSe2 


!wSe2 


1 


1st smuon 


~e2 


-E2 


1000013 





MSmul 


!wSmul 


1 


2nd smuon 


~e5 


-E5 


2000013 





MSmu2 


!wSmu2 


1 


1st stau 


~e3 


-E3 


1000015 





MStaul 


IwStaul 


1 


2nd stau 


~e6 


-E6 


2000015 





MStau2 


!wStau2 


1 


e-sneutrino 


-nl 


-Nl 


1000012 





MSnl 


IwSnl 


1 


mu-sneutrino 


-n2 


-N2 


1000014 





MSn2 


!wSn2 


1 


tau-sneutrino 


~n3 


-N3 


1000016 





MSn3 


!wSn3 


1 


1st u-squark 


-ul 


-Ul 


1000002 





MSul 


IwSul 


3 


2nd u-squark 


~u2 


-U2 


2000002 





MSu2 


!wSu2 


3 


1st d-squark 


-dl 


-Dl 


1000001 





MSdl 


IwSdl 


3 


2nd d-squark 


~d2 


-D2 


2000001 





MSd2 


!wSd2 


3 


1st c-squark 


-cl 


-Cl 


1000004 





MScl 


IwScl 


3 


2nd c-squark 


~c2 


-C2 


2000004 





MSc2 


!wSc2 


3 


1st s-squark 


-si 


-SI 


1000003 





MSsl 


IwSsl 


3 


2nd s-squark 


~s2 


-S2 


2000003 





MSs2 


!wSs2 


3 


1st t-squark 


~tl 


-Tl 


1000006 





MStl 


IwStl 


3 


2nd t-squark 


~t2 


-T2 


2000006 





MSt2 


!wSt2 


3 


1st b-squark 


~bl 


-Bl 


1000005 





MSbl 


IwSbl 


3 


2nd b-squark 


~b2 


-B2 


2000005 





MSb2 


!wSb2 


3 



TABLE VII. Particle content in the E6SSM- 



12.02 



with CalcHEP 



naming conventions and properties. 
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Name 


Value 


Comment 


gib 


0.073 


U(l)-primed coupling with sqrt(l/40) 


hL 


0.393 


Yukawa coupling for S Hu Hd 


hL22 


-0.000357 Yukawa coupling for S Hd2 Hu2 


hL21 


0.04 


Yukawa coupling for S Hd2 Hul 


hL12 


0.0321 


Yukawa coupling for S Hdl Hu2 


hLll 


0.000714 


Yukawa coupling for S Hdl Hul 


hFd22 


0.001 


Yukawa coupling for S2 Hd Hu2 


hFd21 


0.6844 


Yukawa coupling for S2 Hd Hul 


hFdl2 


0.65 


Yukawa coupling for SI Hd Hu2 


hFdll 


0.001 


Yukawa coupling for SI Hd Hul 


hFu22 


0.001 


Yukawa coupling for S2 Hu Hd2 


hFu21 


0.67 


Yukawa coupling for S2 Hu Hdl 


hFul2 


0.64 


Yukawa coupling for SI Hu Hd2 


hFull 


0.001 


Yukawa coupling for SI Hu Hdl 


hXd2 


0.000714 


Yukawa coupling for S Hd Hu2 


hXdl 


0.000714 


Yukawa coupling for S Hd Hul 


hXu2 


0.000714 


Yukawa coupling for S Hd2 Hu 


hXul 


0.000714 


Yukawa coupling for S Hdl Hu 


hZ2 


0.001 


Yukawa coupling for S2 Hd Hu 


hZl 


0.001 


Yukawa coupling for SI Hd Hu 


Svev 


5180 


SM-singlet VEV 


Hvev 


246 


SM Higgs VEV 


MSq 


2000 


Common soft squark mass scale for gen. 1 and 2 = Mql= Mq2= Mul= Mu2= Mdl= Md2 


topA 


-2200 


soft trilinear A-coupling for top 


botA 


-2200 


soft trilinear A-coupling for bottom 


Mq3 


2000 


Soft squark mass for third gen. SU(2) doublet, q 


Mu3 


2000 


Soft squark mass for third gen. right-handed u 


Md3 


2000 


Soft squark mass for third gen. right-handed d 


Mil 


2000 


Soft slepton mass for 1st gen. SU(2) doublet, L 


M12 


2000 


Soft slepton mass for 2nd gen. SU(2) doublet, L 


M13 


2000 


Soft slepton mass for 3rd gen. SU(2) doublet, L 


Mrl 


2000 


Soft slepton mass for 1st gen. right-handed e (selectron) 


Mr2 


2000 


Soft slepton mass for 2nd gen. right-handed e (smuon) 


Mr3 


2040 


Soft slepton mass for 3rd gen. right-handed e (stau) 


tauA 


-2200 


soft trilinear A-coupling for tau 


lsc 


165 


scale for Higgs loop corrections 


lmt 


165 


top mass in Higgs loop corrections 


Mha 


2736 


pseudoscalar Higgs mass 


MG1 


150 


Soft gaugino mass for U(l) (hypercharge) 


MG2 


300 


Soft gaugino mass for SU(2) 


MGlb 


151 


Soft gaugino mass for U(l)' 


Maux 


1 


mass of aux particles 


tb 


1.5 


tangent beta 


MGo 


800 


gluino mass 



TABLE VIII. Input parameters for the E6SSM-12.02 model in CalcHEP notation. Some SM parameters have been removed 
from this list. This is the format and content of the varsNN.mdl file that comes with the model. By default the parameter 
values of benchmark EeSSM-I are given. 

Name Expression 



MNE13 -MZ*SW*cb 
MNE14 MZ*SW*sb 



MNol MassArray(NeDiag, 1) % Neutralino mass 1 
MNo2 MassArray(NeDiag, 2) % Neutralino mass 2 



TABLE IX. Example lines from the model file funcNN.mdl, where dependent parameters are specified. Comments are allowed 
at the ends of lines after a %. The lines shown are examples of simple expressions for matrix elements and masses defined by 
external numerical routines. 



35 



PI P2 P3 P4 Factor dLagrangian/ dA(pl) dA(p2) dA(p3) 

A H+ H- T5E ml.p2-ml.p3 

G W+ -CI -bl -EE*GG*Sqrt2*Vcb*Zd33/SW ml.m2 

W+ W+ W- W- -EE~2/SW~2 ml.m4*m2.m3-2*ml.m2*m3. m4+ml.m3*m2.m4 

TABLE X. Example lines from the model file lgrngNN.mdl, where Feynman rules for all vertices in the model are listed. 



4. Feynman rules 



All of the vertices in the model are listed in the file lgrngNN.mdl. Some of the vertices include auxiliary particles, 
which are included for technical reasons, e.g. to construct a four-gluon vertex. Examples of vertices and Feynman 
rules from this file are given in Tab. [X] 
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